Representative Pathways of Excitation Migration in Photosystem I 
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Photosystem I is a protein-pigment complex that performs photosynthesis in plants, green algae, 
and cyanobacteria. It contains an aggregate of chlorophylls that absorbs light and delivers the 
resulting electronic excitation to the special pair of chlorophylls where the excitation energy is used 
for producing charge separation across the cell membrane. The seemingly random arrangement 
of chlorophylls in photosystem I poses the question which pathways excitation migration follows 
towards the special pair after absorption of light at any of its chlorophylls. We employ a master 
equation to describe the process of excitation migration in photosystem I, and find representative 
paths of excitation migration based on the mean first-passage time from each chlorophyll to the 
special pair. The resulting paths, beginning at each of the chlorophylls and ending at the special 
pair, provide a complete yet distilled picture of the excitation migration towards the special pair. 
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Introduction 

Photosynthesis is carried out by pigment-protein com- 
plexes embedded in cell membranes. An aggregate of 
interacting pigments held in a fixed arrangement in such 
a complex absorbs light and delivers the resulting elec- 
tronic excitation to a reaction center, which uses the ex- 
citation energy to produce a charge separation across the 
cell membrane . The transmembrane potential induced 
by this charge separation is later used for the synthesis 
of ATP. 

Among the photosynthetic complexes, the photosyn- 
thetic unit of purple bacteria has been most exten- 
sively studied; for a review see [Q. Recently, a high- 
resolution structure of photosystem I (PSI) has been ob- 
tained from cyanobacterium Synechococcus elongatus . 
PSI along with photosystem II constitutes the machin- 
ery of oxygenic photosynthesis in plants, green algae, and 
cyanobacteria. PSI contains an aggregate of 96 chloro- 
phylls, including the special pair P700, where charge sep- 
aration is initiated. The chlorophylls are arranged with- 
out an apparent order, except for an apparent pseudo- 
C2 symmetry, which is in contrast to the highly sym- 
metric arrangement of bacteriochlorophylls in the pho- 
tosynthetic unit of purple bacteria. The rather ran- 
dom arrangement of chlorophylls in PSI poses the ques- 
tion which pathways excitation migration follows towards 
P700 after any of its chlorophylls has absorbed light. 

The rates of inter-chlorophyll excitation transfer in PSI 
have been calculated based on Forster theory 0, ||, 1). 
Excitation follows a stochastic trajectory along the exci- 
tation transfer network given by these rates. However, 
the obscure pattern of the excitation transfer network in 
PSI (cf. Fig. 6 in Ref. instills a need for a simpler 
and more distilled picture of excitation migration. Also, 
the aspect that the special pair P700 is the target of the 
excitation migration process is not incorporated in the 
excitation transfer network solely given by the transfer 
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rates. Therefore, we seek paths directed towards P700 
that are representative of all events of excitation migra- 
tion towards P700 in a way that actual migration tra- 
jectories can be considered as noisy trajectories around 
those representative paths. 

The task of finding such representative pathways is also 
encountered in the so-called reaction-path problem |Q, |j . 
Recently, a new approach to the reaction-path problem 
was suggested [||, where mean first- passage times were 
employed as reaction coordinates and reaction paths were 
constructed from them. The suggested method is appli- 
cable to discrete systems such as excitation migration in 
photosynthetic complexes. In this letter we apply the 
method to find excitation migration pathways in PSI. 



Theory 

Excitation migration in PSI is a stochastic process gov- 
erned by the rates of inter-chlorophyll transfer, dissipa- 
tion (internal conversion of excitation to heat), and the 
charge separation at P700. We build a master equation 
describing the process, calculate the mean first-passage 
time for excitation migration from each chlorophyll to 
P700, and then find pathways representative of the ex- 
citation migration in the excitation transfer network of 
the chlorophyll aggregate towards the P700 pair. 

The process of excitation migration can be described 
in terms of the probability Pij(t) that an excitation, ini- 
tiated by light absorption at chlorophyll j, resides at 
chlorophyll i after time interval t. The inter-chlorophyll 
transfer rates Ti^j from chlorophyll i to chlorophyll j 
are calculated as explained in Ref. Q. For the calcula- 
tion of Ti^j , we use the recently obtained site energies for 
chlorophylls ^ and the inter-chlorophyll electronic cou- 
plings given by a full Coulomb computation that includes 
all multipole contributions to the coupling Q . The dissi- 
pation rate is assumed to be the same (/cdiss ~ (1 ns)^^) at 
all the chlorophylls. Since we are here interested in paths 
along which excitation migrates towards P700, it is sufh- 
cient to consider only first passages of excitation to P700. 
(The first passage corresponds to the first term in the so- 
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called sojourn expansion whose higher-order terms cover 
subsequent escapes and returns to P700 [^.) Thus, the 
charge-separation rate is not needed in the model, and 
the transfer rates from P700 to the other chlorophylls 
can be set to zero. Collecting these rates, we obtain the 
master equation ||l0[] 



dt 



(1) 



where the 96x96 transition matrix A';, is 
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From the master equation (Eq. yj) we derive a backward 
master equation which is used for calculating the mean 
first-passage time [0 |l2| . Since excitation migration is a 
Markov process, the probability (t) can be expressed 
as 



(3) 



where = X^iePTOo ^/3-»i the transfer rate to P700 
from chlorophyll p. Thus, the (conditional) mean first- 
passage time Ta , namely the average exit time given that 
the exit has occurred through a delivery to P700 (not 
through dissipation), is given as 
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Here (jja is the total probability that the excitation is 
eventually delivered to P700. Since the excitation may be 
dissipated before it reaches P700, (pa is less than one and 
the denominator in Eq. ^ is necessary for normalization. 
From Eqs. ^ ||, and H, we obtain equations for (ji^ and 
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for an arbitrary intermediate time s (0 < s < t). 
the derivative with respect to t leads to 
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We now take the limit of s 
backward master equation 



(t — 0) and obtain the 



dt 
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where we have used lim^^ 
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P^l(t 



In the model constructed above, excitation can exit 
the system by two means: the delivery to P700 and dis- 
sipation. We emphasize again that we only consider first 
passages to P700. Thus, it is convenient to consider a 
subsystem of 94 chlorophylls, excluding the P700 pair of 
chlorophylls. The subsystem is described by the master 
equation 



dt ^7 
or the backward master equation 



dt 
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where Greek subscripts indicate that the P700 chloro- 
phylls are not included. The 94 x 94 transition matrix 
Kaf3 is obtained from the original 96 x 96 transition ma- 
trix Kij by eliminating the rows and columns belonging 
to the P700 chlorophylls. 

Consider an excitation initiated at chlorophyll a at 
time zero. The probability that the excitation reaches 
P700 between time t and t+dt is equal to J^p ^pPpa (t), 
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The solutions to these equations are given in terms of the 
inverse matrix of the matrix Kap'- 
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The mean first-passage time (tj =0 for i GP700) is a 
natural measure of the progress in the excitation migra- 
tion towards the P700 pair chlorophylls associated 
with shorter mean first-passage times are closer to the 
target (P700) in the migration process. As suggested in 
Ref. , we construct excitation migration paths accord- 
ing to the scheme that a path going through chlorophyll 
i chooses the next chlorophyll j such that Ti^j{Ti—Tj) is 
maximized. The resulting paths can be considered rep- 
resentative because the mean first-passage time Ti is an 
average over all events of the excitation migration from 
chlorophyll i to P700. For a transfer step from chloro- 
phyll i to chlorophyll the quantity 1/Ti^j may be in- 
terpreted as a cost, and Ti — Tj as a gain. The scheme 
then amounts to maximizing the gain-cost ratio for each 
step. 



Results and Discussion 

The above method yields the mean first-passage 
time and the paths shown in Fig. |l|. The 96 chloro- 
phylls are sorted in the order of increasing mean 
first-passage time and are listed in Fig. p](c). The 
excitation migration paths shown in Fig. Wa) exhibits 
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FIG. 1: Mean first-passage time and pathways of excitation migration in PSI. (a) Tfie diloropliylls, projected onto tfie 
membrane plane, are denoted by circles color-coded according to the mean first-passage time to the P700 pair. The excitation 
migration paths constructed based on the mean first-passage time are shown as arrows, (b) The mean first-passage time is 
plotted in the increasing order. The color-code scheme is the same as in (a), (c) The chlorophylls are sorted in the order of 
increasing mean first-passage time. The nomenclature (eC-Al, etc) follows Rcf. [^. 'Red chlorophyll' candidates are marked 
with red dots in (a) and (c). 
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a network without an apparent order, as expected 
from the arrangement of the chlorophylls. There are 
various paths, including six direct paths to P700 (from 
chlorophylls eC-A2, eC-B2, A24, A26, B5, and B24) 
and the most complicated path composed of nine 
steps (A21^PLl^A22^A30^B37^B38^B39->eC- 
B3^eC-A2^eC-Bl). We here emphasize again that 
excitation does not necessarily follow these paths, but 
that they should be understood as representative paths. 

As can be seen from Fig. |](b), most of the chloro- 
phylls (9th to 96th) have mean first-passage times around 
or above 10 ps, with the average over all being 12.5 ps. 
Chlorophylls A40 and B39 lie between these peripheral 
chlorophylls and the reaction-center chlorophylls which 
coincide with the six chlorophylls with the shortest mean 
first-passage times. This supports to a certain extent the 
suggestion that chlorophylls A40 and B39 connect the 
reaction-center chlorophylls and the peripheral chloro- 
phylls 1^. But, among the found excitation migration 
paths to P700, many go through neither chlorophyll A40 
nor B39. Therefore, these two chlorophylls are not the 
only connection between the reaction center and the pe- 
riphery, and should not be considered as bottlenecks. 

The last chlorophyll in Fig. ^(c), chlorophyll Ml, raises 
a question as its mean first-passage time (18.8 ps) is much 
longer than the second longest one (15.6 ps). This docs 
not imply that excitation is trapped at this chlorophyll; 
18.8 ps is still much shorter than the dissipation time of 
Ins. But, it seems inappropriate that one chlorophyll is 
located relatively farther from the rest. PSI exists as a 
trimer in vivo |l^ and chlorophyll Ml lies close to the 
boundary between monomers. By computing the elec- 
tronic couplings between chlorophyll Ml and the chloro- 
phylls in the next monomer, we find that chlorophyll Ml 
is coupled to chlorophyll A30 in the next monomer with 



the coupling of 52.9 cm~^, which is much stronger than 
the strongest coupling it has within its own monomer 
(6.7 cm^^ with B8). Chlorophyll Ml functionally belongs 
to the next monomer, not its own monomer. 

In PSI, it is experimentally known that seven to eleven 
chlorophylls absorb light at longer wavelengths than the 
special pair P700 |lj, These 'red chlorophylls' have 
not been identified with certainty, but some candidates 
(marked with red dots in Figs. |l|(a) and (c)) have been 
proposed ^ ||, |^. Most of the candidates arc found 
to be located rather close to P700 in terms of not only 
spatial distance but also the mean first-passage time. 
Chlorophylls A38, A39, B37, and B38 have the shortest 
mean first-passage times among peripheral chlorophylls 
(except for B39 and A40). With regard to the represen- 
tative paths shown in Fig. |^(a), chlorophylls A24, A32, 
A38, B22, B24, and B26 are one step away from the 
reaction-center chlorophylls, and chlorophylls A39, B6, 
and B38 are one step away from the 'connecting chloro- 
phylls', A40 and B39. 

In conclusion, we have found representative paths of 
excitation migration in PSI based on mean first-passage 
times. The mean first-passage time and the paths provide 
a complete yet distilled picture of the excitation migra- 
tion towards P700 as illustrated in Fig. |l|. We expect 
that our methodology will be useful for various light- 
harvesting complexes as more of their high-resolution 
structures become available. 
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